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The Kondo lattice model, augmented by a Zeeman term, serves as a useful model of a Kondo 
insulator in an applied magnetic field. A variational mean field analysis of this system on a square 
lattice, backed up by quantum Monte Carlo calculations, reveals an interesting separation of mag- 
netic field scales. For Zeeman energy comparable to the Kondo energy, the spin gap closes and the 
system develops transverse staggered magnetic order. The charge gap, however, survives up to a 
higher hybridization energy scale, at which point the canted antiferromagnetism is suppressed and 
the system becomes metallic. Quantum Monte Carlo simulations support this mean field scenario. 
An interesting rearrangement of spectral weight with magnetic field is found. 



The unusual properties of heavy fermion materials are 
the result of strong interactions between their mobile 
spd- and tightly bound /-electrons. In the standard pic- 
ture 1], a broad conduction band is intersected by a 
nearly flat band of renormalized core levels. Hybridiza- 
tion produces a band structure with extremely shallow 
dispersion near the band edge. When the chemical po- 
tential lies below the hybridization gap, the system is 
a metal characterized by a very large effective electron 
mass. When the chemical potential lies inside the gap, 
the system is an insulator. Since RKKY- mediated anti- 
ferromagnetism competes with this hybridization mech- 
anism, heavy fermion systems typically live at the edge 
of magnetic instability [2|. 

The Kondo lattice model (KLM) is thought to provide 
an approximate description of the heavy fermion mate- 
rials. In particular, at half-filling, local singlet forma- 
tion dominates whenever the Kondo coupling J exceeds 
some critical value J c ; RKKY antiferromagnetism wins 
out otherwise. The large J regime is believed to describe 
such Kondo insulators as CesBi4Pt3, CeRhAs, CeRhSb, 
YbBi 2 , and SmB 6 . 

An applied magnetic field will interfere with the 
singlet-RKKY competition and, if strong enough, tilt the 
balance in favor of the antiferromagnetism. The Zeeman 
term, a field B coupled to the total magnetic moment, 
favors triplet rather than singlet formation at each site. 
Accordingly, it suppresses the singlet amplitude and thus 
has the potential to stabilize an antiferromagnetic phase, 
even in the J > J c region. Moreover, the Zeeman split- 
ting lifts the degeneracy of the spin up and spin down 
bands, shifting them with respect to one another and po- 
tentially closing the charge gap. Indeed, transport mea- 
surements of Ce3Bi4Pt3 in high magnetic fields 0, Q 
indicate that its resistivity plummets to a low impurity- 
dominated value at a critical field on the order of 50 T. 

In this Letter, we show by variational mean field and 
quantum Monte Carlo (QMC) calculations that as the 
applied field is ramped up the ground state of the KLM 



with Zeeman splitting (KLM+Z) transforms from Kondo 
insulator to conventional metal by way of a canted anti- 
ferromagnetic phase. 

The KLM+Z describes a half-filled conduction band 
interacting with a lattice of localized spins {Si}. The net 
magnetic moment at each site is coupled to an external 
field. The Hamiltonian is H = H t + HjB, where 



(ij) 



(i) 



Here, t, J, and B are the hopping, Kondo coupling, and 
magnetic field, respectively, (c) is the creation (an- 
nihilation) operator for the conduction electrons. The 
lattice is taken to be square — its bipartite nature being 
essential to our QMC scheme. 

Variational Mean Field — Take S = \pcrf , subject to 
the constraint pf = l : where p (/) is the creation (an- 
nihilation) operator of a fictitious, dispersionless /-band. 
At the mean field level, it is sufficient to enforce the con- 
straint on average and require only that (Pf) = 1. In the 
/ language, the local spin singlet operator is pc. 

A comprehensive variational calculation of the ground 
state must incorporate all of the following: hybridiza- 
tion between the c and / bands (singlet formation), an- 
tiferromagnetism transverse to the field, and magnetism 
parallel to the field. A suitable trial wavefunction is the 
ground state of the mean field Hamiltonian 

Hmf = H t ~Y,{ V *fi c i + Vc Ui) 

i 
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It contains five symmetry-breaking terms, controlled by 
the variational parameters {V m } = {V, M c . Mf, B Cl Bf}. 



2 



B/t 



-1 



-2 



Canted 
b C 2/ AFM 




Metal 



Hybridization 
-1.0 

0.8 



0.6 
0.4 
0.2 



[JJ] [00] [10] 



_ 



Charge Gap 
n >0.5 

0.4 
"0.3 

0.2 
0.1 




_ 



Uniform Mag (c) 

n0.5 
-°oi + 

0.2 
0.1 




Staggered Mag (c) 
— 0.15 
lo.12 

0.09 
0.06 
0.03 




2 

J/t 



FIG. 1: T = mean field, square lattice. (Left) The sin- 
glet phase has hybridization order only. Antiferromagnetic 
and hybridization order coexist in the striped region. The 
stippling indicates that the electron moments are directed 
opposite to the B field. In the metallic phase, the local 
moments are saturated. (Right) (f^c), A c , (^c\<j 3 Ci), and 
{( — ly^cl^a) in the upper half of the phase diagram. 



The variational ground state energy W[{V^}] is the ex- 
pectation value of the exact Hamiltonian evaluated in the 
ground state of Hmf- {V m } is then chosen such that U is 
minimized. Figure^depicts the resulting phase diagram. 

The singlet regime is characterized by nonzero hy- 
bridization (V 7^ 0) and the absence of magnetic order 
(M c j = 0, B c j = B). The energy levels, given by 
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with Eh — — 2tJ2 l=1 cos ki (for a d-hypercubic lattice), 
are parameterized by the band index n = ± and spin pro- 
jection s = ± (T,!)- In the B = case (inset Fig. EJ, the 
band separation takes its minimum, 2V, on the surface 
| fci| + • • • + \kd\ = 7r (ek = 0). The band gap, however, 
is indirect: promoting a quasiparticle from the top of 
the lower band [at Q = (ir, . . . , 7r)] to an arbitrary mo- 
mentum state in the upper band costs uj qp (fc) = E£ —Eq\ 
around its minimum, c^ qp (fc) ^2Ak + (l/2m*)|fc| 2 , where 
41/ 2 = 2A K (2A K +V^) and W = 4dt is the noninteracting 
bandwidth. The so-called Kondo energy, Ak(< V), sets 
the scale for both the charge gap (A c = 2Ak) and the 
ground state energy shift (U[V] — U[0] ~ — Ak). 
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FIG. 2: J/t = 3, B cl /t = 0.323. (Top row) Mean field band 
structure, folded into the reduced Brillouin zone, as a func- 
tion of applied magnetic field. (Middle) Same, magnified to 
show the evolution of the band gap. (Bottom) The greyscale 
(white-black 0-0. 5t) indicates the wavevector-dependent 
gap magnitude. For \B \ <B c i, A c (k) is a minimum at fe = 0, 
and for \B \ >B c i, at an expanding ring of points. 



Now consider Bj^O. Throughout the singlet phase, V 
is independent of the applied field. Thus, according to 
Eq. (jnj), the hybridized bands simply shift with respect 
to one another in response to the applied field and the 
charge gap decreases monotonically: A c = 2Ak— \B\. Be- 
fore the charge gap closes completely, however, magnetic 
order sets in (\B\ = B c i < 2Ak). The localized spins 
develop a uniform moment directed with the field and 
a staggered moment perpendicular to it. The electrons 
do the same but are "effectively diamagnetic," canting 
against the field. As \B\ increases, V begins to falls off 
and vanishes at \B\=B'. When \B\ =B" ^max( J/4, B'), 
(^cJct 3 q) changes sign. When \B\ = J? C 2, the local 
spins become saturated and the antiferromagnetism dies 
out; the conduction electrons are decoupled from the 
quenched local spins and propagate freely. 

Naively, one might have expected the charge gap to 
close at \B\ = 2Ak- As shown in Fig. [2j however, the 
incipient antiferromagnetism prevents the level crossing 
by mixing the (c, f)^ and (c, f)k+Qi bands. The charge 
gap has the behavior shown in Fig. [I] (sideset, right): it 
decreases linearly with \B\ in the singlet phase, resurges 
in the canted antiferromagnetic phase, and finally van- 
ishes in the metallic phase. 

An important feature of the metallic transition is that 
the charge gap does not close at the center of the reduced 
(mod Q) Brillouin zone. As the system is driven through 
the antiferromagnetic phase {B c \ < \B\ < B^), the gap 
migrates from k = out to the zone edge before closing. A 
consequence is that the hybridization energy and not the 
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FIG. 3: Analytically continued data from the one-particle 
electron Green's function for J/t = 1.7, /3t = 14 on an 8 x 8 
lattice. Spectral functions A(k, u>) are plotted for wavevectors 
along high- symmetry lines in the Brillouin zone. (Inset) the 
peaks locations (open circles) are superimposed on the mean 
field band structure [Eq. (|3J] with V chosen to fit the band 
splitting at k = (jr, 0). The greyscale indicates the spectral 
weight (white-black <-> 0-1). 

Kondo energy determines the robustness of the insulating 
state. This leads to a separation of energy scales: the 
spin gap closes when \B \ ~2Ak, whereas the charge gap 
closes when \B\~y/2A K W. 

Quantum Monte Carlo — The partition function of the 
system can be written as 

Z = J drjdfje^irjl JJ P J V Xk e~ e]fl[xk] P\-r)), (4) 

k=l 

where k labels the slices of discretized imaginary time 
(e = /3/L), the Grassman vector rj=(c^ c± /j fi) holds 
the four fermionic species, and the operator P projects 
out states that violate the single occupancy requirement 
of the / electrons. H[xk] is bilinear and related to the 
Hamiltonian by a Hubbard- Stratonovich decomposition 
of the Kondo interaction in the hybridization channel. 

Replacing e~ eH ^- Xk ^ by its coherent state representa- 
tion puts a set of Grassman states at each time slice 
k. The matrix elements (r7fc|.P| 77^+1) can be handled by 
introducing a U(l) gauge field zj^k+i living on the tem- 
poral links. Integrating out the Grassman fields gives 
Z = JV X Vz e~ s ^ , where S = £\ fc \ X i,k | 2 -Tr In M[ X , z] 
is the action of a lattice gauge theory in d+1 dimensions. 



By exploiting the particle-hole symmetry of the Hamilto- 
nian and the bipartite nature of the lattice, it is possible 
to transform to a gauge in which M is positive definite. 
Specifically, under c A -> c iU a± -> (-l)*Ci 2 , M -> fa, 
fil — > (— l) i+1 /i2, the gauge theory acquires a block diag- 
onal form, M diag(Mi, Af 2 ) with M 2 Ml, so that 
detM = |detMi| 2 . Since the total magnetic moment 
transforms as ^ i (c i o- 3 c i + fid 3 fi) ^Y^i Z) s =i,2 (^sQs + 
fas f is), the positivity of M is preserved even for B^O. 

A general n-particle correlation function of the form 
Z~ x \VxVz Mr 1 , , • • . Mr), mrr . e ~ s ^ is evaluated 
via stochastic sampling 5] by interpreting V[x,z] = 
Z~ 1 e~ s \- X,z } as a probability weight. (Since detM > 0, 
there is no fermion sign problem.) The sampling algo- 
rithm ensures that any particular configuration x, z is vis- 
ited with probability V[x,z\. Thus, computing correla- 
tion functions amounts to binning and tabulating Mr}. TT , 
for a large series of independent field configurations. Up- 
dates are effected by a more sophisticated version of the 
algorithm introduced by Blankenbecler, Scalapino, and 
Sugar . Implementation details and other technical is- 
sues are discussed at length elsewhere |jj . Our method is 
related to that of Caponi and Assaad, who have carried 
out QMC for the KLM in zero field Q. The reliability 
of our code is verified by comparison with their results: 
e.g., our computed critical coupling, J c /t = 1.47±0.08, is 
consistent with their value of 1.45±0.05. 

The one-particle electron Green's function G(fe,r) = 
< T [cfc(0)4(r)]> corresponds to fV X Vz M^ T P\x, z]. 
Written in terms of its spectral function A(k,(j), the 
Green's function has the form of a linear functional 
G(fc,r) = JdujK(r,u)A(k,u) = K[A(k,u)] with kernel 
K(T,u) = e- UT /(e- l3 " + l). To extract A(k,u), we per- 
form a (systematically regularized, see Ref. Q) functional 
inversion A(k,u) = K -1 [G(fc, r)]. Figure shows the 
zero-field spectrum in the singlet phase (J > J c ). Note 
that the spectral peaks trace out two bands separated 
by a small gap at the chemical potential. The lower and 
upper bands exhibit (heavy fermion) regions of low spec- 
tral weight and nearly flat dispersion in the vicinity of 
k = (7r, 7r) and k = 0, respectively. 

QMC results confirm that the KLM+Z system is well- 
described by its mean field theory. Several signature 
features are observed in the simulations: a robust sin- 
glet phase (J > J c , I B | < B c \), electronic moments di- 
rected against the field (B c i < \B\ < B"), transverse 
staggered magnetic order {B c \ < \B\ < B C 2), and a high- 
field insulator-to- metal transition (\B\ = B C 2) in which 
the charge gap closes near the reduced zone edge. 

The inset in Fig. 21 depicts the phase boundary be- 
tween the singlet and canted antiferromagnetic phases; 
it also marks where (^cJ<j 3 q) changes sign. The main 
plot shows in six panels the spectral function of a spin 
up electron evaluated at a series of wavevectors snaking 
through the Brillouin zone from k = to k = (7r/2, 7r/2); 
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the panels are arranged (from bottom to top) in order of 
increasing distance from the zone center. 

As per Eq. Q , the spectra exhibit a double peak struc- 
ture and drift leftward (lowering energy) as B increases 
from zero. For small fc, there is a primary peak near the 
non-interacting particle energy uj = Sk and a secondary 
peak near uo = Ak- The leftward drift of the secondary 
peak is interrupted by the growth of antiferromagnetic 
correlations that protect the gap. The spectral weight re- 
arranges by hopping over the forbidden region. There is 
no weight at the chemical potential, so the system re- 
mains insulating even when |£?| = 2Ak. 

At k = (7r/2,7r/2), there are two equally weighted 
peaks at ±V. As B is ramped up, spectral weight from 
the lower peak is transferred to the upper, unoccupied 
peak (this accounts for (^cV 3 c) < 0), which grows in- 
creasingly sharper as it drifts leftward and crosses the 
chemical potential (fixed at \i = by particle-hole sym- 
metry). On the approach to the metallic transition, 
the spectral weight of the lower peak is exhausted and 
A(k, uo) ^5(uj — B/2), its free electron value; the antifer- 
romagnetic correlations die out here as well. Contrary 
to mean field prediction, the closing of the gap does not 
coincide exactly with saturation of the local moments. 

Summary — We have investigated the effect of an ap- 
plied magnetic field on the Kondo insulator ground state, 
using mean field calculations (in the thermodynamic 
limit) and QMC simulations (on small lattices) to charac- 
terize the phases of the KLM+Z. From these numerical 
results, a consistent picture emerges: (i) In the J > J c 
local singlet phase, the hybridized bands respond to the 
applied field by shifting with respect to one another and 
reducing the charge gap. (ii) For \B\ ~ 2Ak, the spin 
gap closes and the local singlet phase gives way to an 
antiferromagnetic phase in which the local spins cant 
with and the electron moments cant against the applied 
field. The hybridization-supported charge gap, which 
had nearly closed, reinflates, driven now by the magnetic 
order, (iii) At larger fields, there is a crossover at which 
the electron moments begin to cant with the field, (iv) 
For | B | ~ ^/2AkW, the antiferromagnetism vanishes and 
the system becomes metallic. 
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FIG. 4: J/t = 1.7, f3t = 14, 8x8 lattice. Spectral functions 
A^(k,u) are plotted over a range of field values (0 < B/t < 
1.8 offset) for a series of wavevectors (identified in the top 
right inset). Smaller peaks have been scaled by the indicated 
magnification factor. Note that, instead of shifting smoothly 
through zero in a magnetic field, the spectral weight jumps 
across the gap, except for k on the reduced zone boundary. 
(Top left inset) A partial phase diagram. Electron moments 
are directed against the field between B c i and B" . 
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